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TECHNICAi:, MEMORANDUM 


LOW DRAG ATTITUDE CONTROL FOR SKYLAB 
ORBITAL LIFETIME EXTENSION 

INTRODUCTION 


On Febi’unry 9, 1974, Skylab systems were confif»'ui’ed for n final 
power down and Skylab was deactivated in a passively stabilissod pa*nvity 
gradient (GO) attitude with the Mutiple Docking Adapter (MDA) up and 
the solar panels trailing (Fig.l). Prediction of solar yycle 21 activity 
(the solar cycle predicted to begin in 1977) indicated that this attitude 
would result in a potential storage period of 8 to 10 yr. However, in 
the fall of 1977 it was determined that Skylab had started to tumble 
rondomly and was experiencing an increased orbital decay rate. This 
was the result of the greater than predicted solar activity at thebeginning 
of solar cycle 21. This inerensed activity incroasod the drag forces on 
the vehicle. Skylab was now predicted to reenter the Earth's atmosphere 
in late 1978 or early 1979 unless i,omothing was done to reduce the drag 
forces acting on it. It was necessary to make a decision to either accept 
an eai’ly uncontrolled reentry (and with it the danger that relatively 
large pieces could damage something or hurt somebody) or to attempt to 
aetively control Skylab in a lower drag attitude thereby extending its 
orbital lifetime until o Space Shuttle mission could effect a boost or 
deorbit maneuver with Skylab, 

In order to verify what options could be accomplished with the on- 
board Skylab systems, in March 1978 a team of NASA engineers went to 
the Bermuda giound station to establish eommunications. The resulting 
data indicated no discernible degi'adation of tlie Skylab systems during 
its four years of orbital storage. The knowledge that Skylab was in an 
unstable tumble prompted investigation into schemes winch might extend 
the orbital lifetime of Skylab. 

The first option investigated was to use the on-board thruster 
attitude control system (TAGS) to maintain a quasistable tumble. How- 
ever, it was soon determined that this option would not extend the life- 
time sufficiently to correspond to the operational readiness of the Space 
Shuttle for a possible reboost or dcorbit mission. The only alternative 
was to I’cactivate and continuously control the Skylab in a minimum drag 
attitude. In order to accomplish this, the low-drag, end-on-velocity- 
vector (EOW) attitude control scheme was developed in record time (the 
authors were p'iven the task on March 20 and the scheme was flown on 
July 11). Skylab remained in the low-drag attitude until January 25, 

1979, when the vehicle was commanded to its original design attitude 
[solar inex'tial (SI)], which was a high-drag attitude (See Skylab EOVV 
Time table for futher details). 




Figui’e 1. Skylab. 


With the active control of Skylab in the low-drag attitude, it was 
decided to accelerate the dovelopaient of an orbital retrieval system (Tele- 
opertor Retrieval System) that might be accommodated on an early flight 
of the Space Shuttle, thus increasing the chances of rendezvousing with 
Skylab. The rate of orbital decay, however, continued to increase due 
to the increased solar activity. Skylab's on-board systems also showed 
signs of deterioration , and there were strong indications that the sche- 
dule of the Space Shuttle would slip. For these reasons, the effort for 
a Skylab I’ecovery was terminated in December 1978, and Skylab was 
placed in the SI attitude in January 1979. In this attitxide the ground 
maintenance was minimum and efforts could be concentrated on a method 
for a conrolled reentry of Skylab [ 1] . 

More detail about the reactivation mission can be found in Refer- 
ence 2. The following section gives the development of the EOVV control 
scheme. 


ATTITUDE AND POINTING CONTROL SYSTEM (APCS) 


The control of the Skylab attitude to the attitude reference was 
done exactly as in the original mission [3]. However, only the pointing 
control system (PCS) of the APCS was used; the experiment pointing 
control system (EPCS) was disabled. 

The major parts of the APCS were the rate gyros, th,* Acquisition 
Sun Sensors (ACQ SS), the Star Tx*acker (it had failed during the 
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SKYIiAB EOVV TIMETABLE 


Authors made aware that a low-draer attitude 
i"omentun mana*yement method was needed. 

Concent review by MSFC middle management. 

Final equations for 2-CMG EOVV to IBM. 

Documentation on ATMDC software change requirements 
for 2-Gf10 EOVV oneration completed by IBM. 

CMG snin up in camed position (7r00 am CST, 
eared to H + CU92 0.03 0.363) 

7 am CDT. Rotation about sunline to place x axis lOF. 

8 am CDT. CMG control in SI attitude, 

12 noon CDT. Loss of control due bo inadverdcnt servo 
power cut-off to CMutfS (faulty switch selector 
introduced additional command) 

1 pm GDT. Rerained control with CMG's after 2 orbits. 

8; 27 dT, Entering EOVV A attitude. 

Loss of EOVV attitude control due, to large angular 
momentum and CMG saturation. 

Re-e$tabilishment of EOVV A attitude,. 

Loss of all Power introduces loss of attitude control. 

Reorientation maneuver to find attitude. 

Re,-establishment of SOW A attitude. 

Maneuver from EOVV A to EOVV F* attitude 

NASA EQ press aonference ( J . Yardley// . Aller) where 
Skvleb is declared unsavaMe. 

Reorientation maneuver to SI attitude (hold mode to 
allow time for the design of the torque equilibrium 
control method for reentry). 
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orifdnol mission) • tho Apollo Telescone Mount Distal Computer (ATMDC), 
the Workshop Computer Interface Unit (WCIU)» three double- gi»*haled 
Control Moment Gyros (CMGs), and cold-pas (compressed nitropen) 
Thruster Attitude Comtrol System (TACS), 

Six control modes were addressable: (1) STANDBY, (2) SOLAR 

INERTIAL (SI), (3) EXPERIMENT POINTING, (4) ATTITUDE HOLD /CMC, 
(S) ATTITUDE HOLD /TAGS, (6) ZLV (for g axis alonp the local vertical), 
EOVV coirtrol was propratamed to be a substate of the ZL’^ mode. The 
basic ZLV attitude was with the positive % axis along the ideal vertical, 
pointing up, and the positive x axis In tho orbital plane, pointing in the 
direction of the velocity vector. Any angular offset form the basic ZLV 
attitude (offset identified by the quaternion Qal) could be commanded by 
a set of three Euler angles (x's) with a y,z,x, rotation sequence. None 
of the original APCS capabilities were eliminated by the addition of the 
EOVV control method. 


CMC CONTROL SYSTEM 


The CMG control system was composed of three orthogonally 
mounted, double gimbnled CMGS with angular momentum magnitude H of 
3050 Nms (2280 ft*lb-soc) as shown in Figure 2. The CMG control law 
utilized three normalized torque commands and the CMC momentun) status 
to generate the proper CMG gimbal rate commands [4]. The ( MG control 
law consisted of three parts; CMG steering law, rotation law, and gimbal 
stop avoidance logic. There also were some other routines for specialized 
situations like caging the CMGs to a desired momentum state [5], 

The CMG control law hod the cap ability to operate with either three 
or two CMG's for redundancy. Since (.MG No.l had failed during the 
original Skylab mission the CMG control law was always in the two-CMG 
option . 


EOVV MOMENTUM CONTROL 


The problem for EOVV was to determine variable reference attitudes 
such that, on the average, the angular momentum was contained within 
the two-CMG capability (allowing the CMGs to hold the presex-ibod attitude 
reference), while the average refei’ence attitude was consistent with the 
desired low aei’odynnmie drag. 

Minimizing tho drag on Skylab requix’ed tJiat the least frontal ax’ea 
was presented to the wind while at tho same time holding a GG toi’que 
equilibrium (at the altitudes of concei’n the GG torques were still very 
dominant and they were therefore used exclusivoly for momentum conti’ol; 
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Figure 2. CMG mounting arrangement, 

for details on the aerodynamic torques, [1]). Keeping the minimum prin- 
cipal moment of inertia axis parallel to the wind direction fulfilled this 
requirement. To have the necessary electrioal power from the solar cells 
as well as stx*apdown update information from the ACQ SS , the Skylab 
was rolled through the angle (APPENDIX F) about the minimum 

principal moraent-of-inertia axis (principal x axis) such that once per 
orbit the sun line passed nominally through the center of the ACQ SS. 

There were two attitudes which satisfied these requirements: One 

with the MDA forward (EOW A) and one with the MDA backward 
(EOVV B). In either case, the MDA had to be pitched down by a varying 
amount (depending on the solar elevation angle 0 “ above the orbital 

plane) to align the principal x axis with the orbit tangent. 
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Control of the ongular momentum was split Into the control of the 
momontum component perpendicular to th*^. orbital plane (POP control) and 
the component in the orbital plane (rOP cL..lrol). Since lOP control hod 
some effect on POP controls lOP control ia ti^atod first. 


lOP Momentum Control 

When Skylab was originally designed, it was desirable to minimize 
the GO torques about the minimum princlpal-moment-of-inertia axis os 
much os possible since the momentum management scheme [63 was least 
efficient about this axis. For EOVV control tins meant that there were 
basically no large GG torques available about this axis and, furthermore, 
there would be no change in torques when the pHnclpul y or z axes were 
±45 degrees from the orbital plane (tracking of tho sun by rolling about 
the X axis would make this a frequent occurrence), Therefore this first 
order effect had to bo abandoned. The actually used momentum control 
scheme assumed that the Skylab momonts-of-inertia were cyclindrical , with 
an average momont"of“inertin difference of AX, 

The problem was solved by using a second order effect. I^rst, a 
large cyclic POP torque was generated by "nodding" (pitching) the Skylab 
in the orbital plane. Cyclic nodding was requii’ed to avoid a continuous 
momentum build-up in POP direction . The cyclic POP torque was the*' 
tilted as required (difisifcntly for each half cycle or quarter orbit) to 
generate a controllable component in the orbital plane. The second order 
effect stems from the fact that the lOP torque is proportional to the nod- 
ding angle times the tilting angle. 

The effectiveness of the lOP momentum control did not depend on 
the frequency of the nodding. However, other considerations entered: 

(1) the lower the frequency, the larger the POP momentum swing, and 

(2) the higher the frequency, the larger the maneuver momentum that has 
to be exchanged between the voWcle and tlic CMG system. Since only a 
limited momentum volume was available with two CMGs and their associated 
ginibal stop iiroblems, a nodding frequency of twice orbital frequency wos 
chosen as a viable compromise (si sin) 
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n s ( 2H t 
''ym 0 




( 1 ) 


Therefore, the lOP momentum control calculations were done every 
quarter orbit. This had the added advantage that the resolution from the 
nearly inertial O system to the rotating L system happened in 90 degree 
intervals allowing indexing of some of the saved momentum samples rather 
than requiring a full-fledged resolution. 

To minimize transients at the quarter orbit sample points, the tiling 
angle was also sinusoidal with twice orbital frequency and its amplitude 
was the only changing quantity; 
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where 
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®a + Oj./4) 




®zo “ n‘omentum component to be dosaturatod 


Ojj - amplitude 


6j^, = ramp per orbit 


The phasing of the nodding and tilting angles was the same so that 
the amplitude was reached half-way between sample points. 

The sample points were chosen so that one of the samples (sample 2) 
occurred at the time when the sun was perpendicular to the solar panels. 
This happened before orbital noon for EOVV A and after noon for EOVV B 
(the difference of about 11 degrees between the geometric and the prin- 
cipal X axes is the reason). In addition the nodding i*ate, being added 
to the ox'bital rate, was phased so that it slowed Skylab down when the 
solar panels were perpendicular to the sun and therefore maximized the 
power from the solar panels [giving rise to the minus signs in equation (1) 
and equation (2)]. 

The tilting angle amplitude was ealculated so that the lOP momentum 
component, which could bo nffcoted during the next quarter orbit (it was 
along the direction of the connecting line between the present sample 
point and the next one), would be driven to the desired value. The 
desired value was basically zero, hut any constant toi’que in the L system 
caused only a cyclic momentum change (with normalized amplitude e^) over 

one orbit and should not be compensated for. Therefore, the momentum 
attributable to a constant E system tox’que was subtracted out of the 
momentum e to be desnturated over the next quarter orbit. 

4uO 
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To recogniiso a cycUo as woU as a I’amp momentum change « four 
past momentum samples were savei. The samples were also used to gen- 
erate stropdown update information once an orbit (flow charts in 
appendix H). 


POP Momentum Control 

The torques ossociated with a rotation about the orbit normal are 
much stronger than the ones associated with lOP control, Hence, the 
momentum sampling for POP control has to bo done as frequently as pos- 
sible. However, the transients should have hod a chance to settle before 
the next POP sample is taken, Twelve somples per orbit satisfied both 
requirements. The desaturation gain (APPENDIX D.b) 


Kyf, ~ H/(3i>^^ A I 


rp 

^des'* 


(3) 


is calculated so that a step attitude change of 


'yc 


Sc 


(■ 1 ) 


desaturates the desired Ae^ in one desaturation interval, To further 

reduce transients, the calculated POP angle (which would have eUminotod 
the momentum offset within the next desaturation interval if the angle were 
applied fully during the interval) was ramped in so that the angle was 
achieved at the end of the interval. Since this only reduced the momentum 
offset by half, the angle was ramped-out during the following interval for 
a full momentum offset elimination. The ramp due to the newly calculated 
POP angle was simpfy superimposed on the ramp-down from the previous 
POP angle (Fig. 3). The ottitnde command is, therefore, given by (N 
signifies the present and N-1 the past value) 


(n 


yc% 




[(ACy)j^ - 0.5 (Aey)j^_^l 


(5) 


The ramps connecting the i] '8 are generated by 
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<"y>N = <"y>N-l + A'V « 


( 0 ) 




Where 

Aly = t(.ly„>N - (llye'N-l’/’-'dos • 


This method resulted in n constant hang-off when necessary: The angle 

change due to the old angle being ramped-out was compensated by the 
ramp-in of the new angle (in fUght, constant hang-offs wore common due 
to strapdown and navigation errors and they were not detrimental, since 
the momentum control kept the voliicle at the truely desired attitude). 

A block diagram of the KOVV orbital y momentum control scheme is shown 
in Figure 4. 


EOVV STRAPDOWN UPDATE 


Strapdown updates about the vehicle x and y axes wore always 
furnished by the ACQ SS. To do that, the roll angle about the principal 
X axis was changed by large amounts to compensate for the large beta 
angle changes (a slow change) and relatively fast smaller corrections were 
applied to compensate for the nodding and the tilting angles. The overall 
effect was that the vehicle z axis nominally traced a cone about orbita’ 
north. The difference between whore the sun-line was at the closest 
approach to the ACQ SS center and where it was supposed to be according 
to the strapdown information gave the strapdown x and y information. 

Since thex’e was no other sensor available, it was more difficult to 
gain strapdown update information about the stin-Une. The selected 
momentum control method, fortunately, had the feature that, due to the 
nodding angle, a misalignment between the ideal orbital plane and the 
indicated orbital plane generated an lOP momentum ramp. 

The actual strapdown update was done by changing the x’efcrence 
quaterion (APPENDIX H.6) 


Qvi = ^ Qyi 


(7) 


where the double-bar operator is defined in APPENDIX E and 


AQ = UQ^, AQg. AQ^l = AQ^l , 


( 8 ) 
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with 


A_q. = 0.5 (s X v + n s) , 

/a 

ond 

AQ^ = /l ~ AQ ■ AQ . 


( 9 ) 


( 10 ) 


The cross product in equation (9) is the ACQ SS update and the 
last term is the lOP ramp update, whore v is a unit vector* in the sun 
direction as calculated by the ATMDC and““ s is the measured sun direction 
unit vector, is the angle about the measured sun direction (the gain 

is developed in APPENDIX D.e.); 

M * 

^°TLN1 “ ‘^TLNIP^ ’ 


where is a gain and ~ ^Xx^NlP) angular momentum 

change (ramp) per orbit (p^ is calculated at sample point 1) modified by 
the ground commanded ramp bias e^,j^ since 


®TLN1P " ^TLNIP ‘^rb ' ^ 

is calculated right after the p calculation and, therefore, is used for 

A 

the next p^ calculation (djij^j^^pis modified at every sample point to 
account for the momentum ehangos commanded by . 


EOVV OPERATION AND PERFORMANCE 


The original EOVV equation considered EOVV A only. In 
EOVV A CMC No. 2 received more solar radiation when the sun was north 
of the orbital plane (positive beta angle) and loss when it was south of 
the orbital plane (negative beta angle). For large negative beta ang:les 
the CMC No. 2 bearing temperatures became critically low. As a con- 
sequence, the EOVV A equations had to be modified duxnng EOVV opei’a- 
tion to allow an EOVV B attitude during extended periods of large nega- 
tive beta angles (Fig. 5). 
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Figure 5. Sun elevation angle beta (April 1978/August 1979). 


Constant angle hang-offs (caused by strapdown, navigation, 
yc 

and other errors) requii’ed constant POP momentum hang-offs to generate 
the necessary commands. Since the range of acceptable POP momentum 
component was rather limited (±0.4H from the nominal; the nominal e^jj 

being 2.5H in EOW A and 0.3H in EOVV B) the nominal POP momentum 
had to be changed to accept large angle hang-offs (3 deg of POP angle 
hang-off required O.IH POP momentum hang-off). This change in nominal 
momentum was made from the ground at the beginning of the EOVV opera- 
tion and later was automated on-board to guard against ground inatten- 
tion, gi'ound system failures, and long telemetry coverage gaps (Fig. 6.). 
Momentum excursions outside the specified range caused loss of attitude 
due to CMC saturation on one occasion June 28, 1978, and it was there- 
fore very important to keep the POP momentum bounded. 

Strapdown updating about the sun- line was done with information 
derived from the lOP momentum ramp. Unfortunately, the evaluation of 
the lOP momentum ramp yielded very noisy readings from one orbit to 
the next and could only be used with a very reduced gain K . This, 
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in turn, could lead to large orbit plane misoligninents to generate the 
required strapdown updates to keep up with the ±50 deg rocking of the 
true orbital plane (due to the precession when viewed with respect to 
the projection of the sun-line into the orbital plane). The large size of 
the maximum change per orbit (1.2 dog) was not I’ocognized at the begin- 
ning of the EOVV operation and was the cause for loss of attitude on 
June 28, 1978. After that the ideal strapdown update necessary to follow 
the rocking of the orbital plane was introduced open loop through the 
quantity called e^,j^ [equation. (12)] and no more problems were experi- 
enced (APPENDIX C). 
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APPENDIX A 


COORDINATE SYSTEMS AND TRANSFORMATIONS 


The coordinate systems which are pertinent to Skylah EOVV control 
are defined here. Each system has some special geometrical or physical 
feature which simplifies the solution of a particular problem. 

The following coordinate systems are described: Principal, Orbital, 

Vehicle, Attitude Reference, Solar Inertial, and Z-Local Vorticol. Each 
coordinate system consists of a set of mutually orthogonal axes exhibiting 
right-handedness . 

An inertial (with respect to rotation only) coordinate system is a 
system which retains its orientation with respect to the celestial sphere, 
although the origin may be moving along any generol curvilinear path in 
space. Similarly, a vehicle fixed system retains its orientation with 
respect to the veliiclo. 


Orbital Coordinate System (O) 

The Orbital Coordinate System (x^, y^, z^) is a processing coor- 
dinate system with its origin at the Earth oenter of mass. The rate of 
precession about the Earth's north polo is approximately -5 degrees /day. 
The z^ axis Ues in the orbital plane, positive through the ascending node 

of the orbit. The x^ axis also lies in the orbital plane 90 degrees ahead 
of the z^ axis. Since the Sky lab orbit was in the x^z^ plane at all times, 
the y^ axis was parallel to the orbital angular momentum vector, complet- 
ing the right-handed system (Fig. A-1). 


Solar Inertial Coordinate System (I) 

The Solar Inertial Coordinate System (Xj, yj, Zj) is only a pseudo- 

inertiul system since it makes one revolution per year. It was used during 
the Skylab mission to point the instruments in the desired direction. The 
oi’igin is coincident with the origin of the Vehicle Coordinate System origin. 
The Zj axis is positive toward the center of the Sun. The Xj axis lies at 

an angle v„ from the orbital plane (this angle is calculated on-board such 
z 

that the principal x axis is in the orbital plane to minimize th? build-up 
of angular momentum) and is positive toward the stsnset terminator. 
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Figure A-1, Coordinate system. 


Z-Local Vertical Coordinate System (L) 

The z-Local Vertical Coordinate System (Xj^, ® rotating 

system with its origin at the center of mass of Skylab (the rote of rota- 
tion is one revolution per orbit). The Xj^ axis is positive in the direction 

of flight and lies in the orbital plane. The Zj^ axis is parallel to the local 

vertical direction and is positive outward, away fi’om the Earth. The yj^ 

axis is parallel to the orbit normal and is positive towai’d orbital North. 


Vehicle Coordinate System (V) 

The Vehicle Coordinate System (Xy, yy, Zy) is a vehicle-fixed 

system with its origin at the center of mass. The Xy axis lies along the 

long axis of Skylab and is positive in the direction of the Multiple Docking 
Adapter (MDA). The Zy axis is positive toward the Apollo Telescope 

Mount (ATM) and the yy axis completes the right-handed system. 
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r^rincipal Axes Coordinate System (P) 

The Principal Axes Coordinate System (Xp, yp, Zp) is a vehicle- 

fixed system v/ith its origin at the center of mass. The axes are along 
the principal moment-of-inertia axes, labeled such that the eigen angle 
between the V and the P system is minimized. 

The following transforms are useful (the subscripts 1, 2, 3 indicate 
rotation about x, y, z, respectively). 


[lil] = [An^ylg 

[AL] -- [Kl^ tnli ir\J^ tHy + %y'l2 
[PV] = [K] 

V 

[VI] = strapdown matrix (updated by w , sunsensor data and 
momentum data) 


APPENDIX B 


GRAVITATIONAL TORQUE MODEL 


The gravitational torques on a satellite produced by a loi’gei spheri- 
col primary body are important contributors to its rotational dynamics. 

The force on a point mass m exerted by the px‘imory M is 


E c - JSKEL R. ; 0 » 6.672E-11 . <BX) 

-* |Rj|3 


The satellite cun be viewed os a collection of point masses. The net tor- 
que on this collection of masses about the origin of satelUte coordinotes 
is 




r. 

HI 



(C2) 


The vector r ^ is the position of m. relative to the origin which is at R 
relative to the primory center. Thus 


R, 


5o*5-i 


(B3) 


and 




r. X GMm. 
—1 1 


R., + r, 



(B4) 


In general |r.j << |R | and hence an expansion of T keeping only low 
order terms becomes sufficient for most purposes. Now 


gp 

|Rp + 




R„ • r. 

3R ^ + h.o.t. 




(B5) 
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Using the expansion of equation (BS) in (B4), 


T 


$ 


r X 
~cm 


GMm 


+ 3 

‘(w»Q 


GM 

iRj' 


n. 


nij rj 


X R, 


R. 


+ h.o.t, . (B6) 

We have used the definition s= njj tj, Rearranging and gi*ouping 
equation (B6) wc obtain 





So + 


3GM 



R. 


o 





+ h.o.t. 


(B7) 


The term in parentheses in equation (B7) occurs in the definition of the 
moment of inertia dyadic (or tensor) 


L = Z (&■ 


1 

ar 


Vt 




(B8) 


Using this definition in equation (B7) and dropping the higher order 
tei'ms j^eld?, the gravity gradient torque expression 


X R + 3 


R, 


•gg ^m |3 ”0 i„ i3 


Is 


-H3 I 


l«ol l«oi 


X I 


A 

l5oi 


(B9) 


As can be seen, if the origin is positioned at the center of mass the more 
familiar gravity gradient torque expression results; 
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I 


(BIO) 


1 


. 3 - 2 “ A. . I . A 

Kr isJ “ IB- 


o' 


For on orbitin^r body m , the orbitol ongulor velocity magnitude 
by 


« GM/(K^ 


Noting that Ro/jR^jl is a unit vector Uj^, we Anally write 


^ J, • Hr 


The torque of equation (B12) is 
gravity gradient torque , 


what is commonly referred to as 


given 

(BID 

(B12) 

the 


APPENDIX C. 


STRAPDOWN DRIFT ERROR AND MOTION OF I 
COORDINATE SYSTEM 


The basic reference coordinate system tised during the original 
Skylab manned mission was the so«called Solar Inertial (1) coordinate 
system (Appendix A gives the definition). As stated there , I is not truly 
an inertial coordinate system since its axes rotate with the sun ond orbit 
regression. Its rotation rate varies with the solar elevation angle out of 
the orbit plane. By definition, the z axis points toward the sun and the 
unit vector Upj^ lies in the orbit plone pointed generally parallel to the 

vehicle velocity vector at orbitol noon. The unit vector Up^ is the 

direction of the x principal axis in veliicle coordinates V. Thus, when 
the V and I systems are aligned, the vehicle x principal axis is in the 
orbit plane so that ideally all gravity gradient torques are cyclic, The 
sun angle Ty and orbit regression angle ,\y are updated once per orbit 

at orbit midnight in the Skylab novigation calculations. In between mid- 
nights, the I system, as defined by veliicle on-board navigation, does 
not rotate. As a consequence of this, it Is useful to redefine I such that 
the definition given previously is only satisfied at midnight. Let Ij^ 

represent the true reference position (we use the convention that true 
or physical parameters are labelled by a subscript k) of this system and 
I represents the estimate. This ostimote is made using the veWcle angular 
rotes as measured by the rate gyros. 


V V V 

“vi “ iJ^VIk ~D 

V V V 

The ongular rates w yj , a P and w are all given in the vehicle 

system V. The rate is tl g>a*o drift rate and experience indicates 

it tends to be constant in the vehicle system V. To a good approximation 
the vehicle remains fixed relative to the local vertical system L. Thus 


V V 

“ VI " % ‘ 


where is the orbital angular rate. 


Since Uj_/ end eonstants, constant 

-Vlk'^ = ^>0 HL 2 '' ■ 2d'^ • 

Wo can integrate the rate 2y,/ using the differential equation 

^ [VIjjl = - Ji (VI|^) . 

This yields 

[VIj^] = exp A [VIj^]^ , 

with 

' n4 ‘ ■ ®‘" 1^1 A + |A|) aA , 

■ |A|^ 

and 

A = - iSvik'' T • 

where 

T = 2n/Q^ , 

is the orbital period. 
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IVt,,]- is the [VI. ] trnnsforimition ut the begimiing of the orbit. 

K O K 

^et 


Svij^ = ^VIj^ T 


and 


then 


2d = 2d ■ 


2vi^ = 2"aL2 ■ 2 d 


iVe shall now assume that |0|j | << 2ir. In practice [0^ | is at most 

V 0.05 to 0.06 radians. Thus wo can evaluato exp (“Syj ^ ] keeping 

. V \ k / 

only 1st order terms m . First, wo find 


I 2ir “* Ur * 0 


-VL 


L2 “D 


Using this result, wo find 


i'VI 


•VI, 


~ -1,2 " Itt “ h} ‘ Hl 2 Hr,2^ 


From this wo find 


V V 

exp (- Oyj ) = exp (Uj ^2 * Id -L2 ^ 
k 


Finally 


tvi^l = exp (Op . Up/) 
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(Cl) 


* (1 + Id ' -L2 Ml2 ^ 


The result we have obtained in equation (Cl) can be suaimarized by say- 
ing the strapdown error induced by gyro drift to the 1 st approximation 
only accumulates along the orbit normal and tends to be purely cyclic 
along Uj^j^ and u^ 3 . 

Strapdown orx’ors enter in also due to motion of the 1 coordinate 
system. These motions are not accounted for in the strapdown calculations 
and so must bo corrected by updates from sun sensor data and momentum 
accumulation data as explained elsewhere. The solar inertial system motion 
is readily calculated fi’om relative motion of the sun and fi*om regression 
of the Skylab orbit. The sun moves approximately 1 leg per day wltilo 
the orbit processes nearly 5.5 deg per day in the retrograde direction. 
Since there are about 16 orbits per day, these effects are usually small 
but occnsicnaliy amount to more than 1 deg per orbit. To analyze these 
effects let us define the system S such that the z axis points to the 
Vernal Equinox and y axis to the eeUptic north pole, 


[SE] = = 


cosfj), sini 0 

?4 Z 

sinii*,, cos* 0 
z z 

0 01 


; ({)„ e: 23.45° 
z 


Also, define R such that Z|^ points to the sun and yj^ to ecliptic north. 
From R to I is a rotation about Zj^ 


[IR] = 


■* [lEl = [O^lj [I'yJj tO.lj , 


where the numeric subscripts indicate the axis of the Euler angle rotation 
( 1 , 2, S represent x, y, z respectively). 
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The angular velocity of I relative to E ie 


w 



0 


• ■ 

sin 6 


• 

z 

0 


cos 


y 

z 

1 

0 


The trons formation [IE] can also be computed from onboard navigation 
parametors 


[IE] = [V 1 [nj [%] [XJ [XI 
Z3 Yg Z3 Y2 


(C2) 


The angle X is the orbit regression angle and X is the nearly constant 

y y 

I'cgression rate. The angles in equation (C2) can all be expressed in 
terms of X^ and The vector Upj^ is a row from the vehicle-to- 

principal axes transformation [PV], 


u 



[PV] 


1 

vH) 

i — 



0 

If 

'■''12 

0 


>■''13 


For convenience, we let [K] = [PV] so that 




11 

12 

13 


The unit vector to the sun is u 


R3 


and is 



sin r coscj) 
y " 

sin r„ sint})„ 

y z 

cos Ty 
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Similarly, the orbit normal is 


-sin A cos A 

« y 

cosA^ 

it 

sinAjj sin Ay 


The angle is by definition the position of orbital noon, i.e., the posi- 

tion in orbit where the radius vector, the solar vector, and the orbit 
normal are coplanar with • Uj^g > 0. 



Hy = tan 


- 1 / ^R3x\ , 

“r3z ' 


whex*e u®* 

IV li 


= fOE] Uj^g^ 


[OE] = tA J [A 1 


riy = tan 


■( 


cosA„ cosA„ sinT., cos* + sinA„ sinP sincj) - cosA„ sinA cosT 
z y y ^ ^ y z it y y 

sinAy sinTy cos^^ cosAy cosFy 


) 


Note: The tan ^ technique used must be a 4 quadrant technique! The 

angle is the elevation angle of the sun out of the orbit plane with 

A 

orbit south the positive hemisphere. 


”x " “R3y 


= sin [sin^^ ®”^‘ y ' ®‘"^z ®*”^y 


The remaining angle v is obtained by solving 

2< 


—PI ’ -02 ~ ^ ** —PI * -LI - ® orbit noon 


(C3) 
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(4 Quadrant tan"^) 

With these intermediate angles and the values of X^,, Ty, (j)^. we 
can compute [IE] from equation (C2) for tony time. We can also compute 

-cosTy cos({)jj cosXy sinXjj + cosTy cosXjj sin<J>^ - sinTy sinXy sinX^ 
sincj) cosX sinX + cosX„ cos({)_ 

i* J it U Z 

-’Sinr cosA GosX sinX + sinT cosX„ sincf) + cosT,, sinX„ sinX 
y ^ y ^ y u 7t yy 55 

From equation (C4) then vvo obtain 


- 

— Oz 


o . -1 13 02z 

f>,,, = sin - 

/(l - (l - Uo«") 

We can now use equation (Cl) to determine Wjg • Figure C-1 shows how 
the components of vary with time as the sun elevation angle (the 
so-called 3 angle 3 = - goes through its maximum value in the north- 
ern orbital hemisphere. We had not realized prior to this that the solar 
inertial frame could move so fast and were unpleasantly surprised when 
EOVV attitude was lost due to strapdown drift and I'esultant momentum 
saturation. For more on this incident see Reference 2. After this inci- 
dent, we developed a procedure to update the z momentum bias in such 
a way that the necessary rotations would be supplied by the bias value 
rather than a momentum error. The parameter was called a ramp bias, 

®rb (12)]. Since the z update angle depends on e^ - e^^, we 

could keep e^, small by using e^^ and since we now knew what to expect 

we could plan in advance what e^^ updates were needed to keep e^, small. 

Table C-1 and Figure C-1 show the e^^ updates that should have been 

made to prevent the loss of momentum control we experienced. The 
knowledge and experience gained here allowed us to successfully pass 
thx’ough two similar peaks of opposite signs in November and January. 
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tan 


■1 ^11 ^02x ^12 ^Ozy 

^11 “02y " ^12 '^02x 




DEG/ORBIT 


1.6 
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TABLE C-1. SCHEDULE OF ERB UPDATE TIMES 


a 




leit 9189 

“23 



161(13144 
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2 

181(81(66 

20 

170t 0113 

1 

188( 2( 9 
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.0 
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‘18 

173116113 

1 

188(10(49 


174113127 

I2 

182(15(32 
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3 

182(20(23 

15 

17SU5I23 


183( K26 

“14 

1761 0134 

j 
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“13 
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6 
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“12 

1761141 3 

"7 

183(19(20 

11 

176119121 

’8 
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“10 
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"9 
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“9 
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184(18(50 

"8 
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“11 
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“7 

177U2I29 

“12 

185(161 9 

“6 

177115141 

“13 

186( 5(37 

*5 

177118152 

“14 

186(21(42 

“4 

I77»22i 4 

“15 

187(17(42 

“3 

1781 11 5 

“16 

188(18(53 

"2 

1781 3t51 

"17 

190( 3(27 

1 

178t 6136 

“18 

191(21(33 

0 

1781 9>22 

19 

194( 0(48 

1 

1781121 6 

20 



178H4I55 

“21 



178«17«43 
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ASSUMES KMUZ IS 

0.2 

178120130 

“23 


178123118 

"24 



179< 2147 

“25 



1791 6130 

26 



I79il0a3 

“27 



179116136 

‘28 



ie0il3tl8 

"27 



180I19U1 

"26 



1811 0147 

"25 



181i 51 8 

“24 
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APPENDIX D. GAIN CALCULATIONS 
a. Calculation of lOP Desaturation Gain K 

nz 

The locol vertical in principal axes components is (s ~ sin, c =: cos) 


cz 

sz 

0 

cy 

0 

-sy 

0 


— m 

-czsy 

-sz 

cz 

0 

0 

1 

0 

0 

= 

szsy 

0 

0 

1 

•sy 

0 

cy 

1 


cy 


where 

y = 'Vn = ' = ’’zn = ' . 

and the rotation is neglected (a cylindrical inertia distribution is 
assumed with I„ = I_ = I and AI = I - I ). The gravity gradient torque 

J *4 /Ik 

is (in P components) 


T 

-ggP 


K 


gg 


■ 0 " 


“'•l ^3 
^’2 

s K 

gg 


0 

czsycy 

-czsyszsy 


where 


- say I 
gg 0 


and in 0 components (where Y = fi^t + y) 


T “ K 
ggo gg 


cY 

0 

sY 

0 

1 

0 

-sY 

0 

cY 


cz 

sz 

0 


-sz 

cz 

0 


0 

0 

1 


0 

czsycy 

"szczsysy 
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K szczsy 


'‘■cy 1 

cy/tz I 

J 


For the further development only the lOP component is of interest, winch 
is (assuming y and z small) 


Sggz.jc = “gg lym "w" 


Sfi^t 


-cJLt 


Integrating over a quarter orbit yields 


H 


■ggoz,x 


ir/2S1^ 

^ ''zni J (s Sl^tc ll^t) 

0 


SSot 




^ ^yni 


ir/2ft 


o 


/ dt sfi t (C^n t - c^^O t) 


ir/2il^ 

“ J dt cVi^t (s^S^^t ~ 
0 


Substitution yields 


(u - cH^t; V = sQ^t) , 


H K K n r> 

Sggoz.x 'gg 'ym "zm 


/du (u^ - u^) 


0 


■f dv (v^ - v^) 
0 


_ 4 


'' ^^yni T5 


1 

-1 
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1 


SggozjX “ TF *^201 


Rotation to the effective axis system yields 


H 


24 , 

ggz,x - IS V* V 


1 

/I 


1 1 
•I 1 


1 

•1 


“ggz " ® ' 


" (24 /T/i5) iMI n T] 
ggx I ' 'ym 


Resolving actual momentum to bo desaturated into the same axes 


1 

’ "z + «x' 

„ i 

1 1 


' »z' 

/I 

-»z * "x 

•» 

/2 

-1 1 


«x 


Only the x component can be desaturated leading to (in units of the 
nominal angular momentum H of one CMC): 


(H//T) (0^ - e^) = (48/15/2) sl^AI Hy„, n,„, 


or 


n 


zm 





> 


with the gain being 

= (15H)/(48!)^AlAy^) 
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b . Calculation of POP Gain K 

yo 

The gravity gi-adient orbital y torque for a i-otation n about the 
orbital y axis through riy is (APPENDIX B> 


^ggoy ” *^gg ®^y 


If the angle iiy is held constant over the desaturation period (about 

1/12 of an orbit) wo get an orbital y momentum change of 

^®ggoy " ^'“gg '■‘■’des ’ 

in units of the nominal CMG momentum H. For a given Ac we need 

oy 


with 


c. Calculation of Strapdown Update Gain 

The components of tlie local vertical in the P system are (neglecting 
the tilt angle n • $ = sin, c = cos) 

Zlli 


cN 

0 

-sn" 


1 

z 

o" 

cy 

0 

sy 
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-sn 
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1 
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~z 
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0 

0 
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0 

= 

-zsy 

SN 

0 
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0 

1 

~sy 

0 

cy 


1 
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with 


N « y + nj y » ft^t; ii « - ny„jS(2il^t){ sPj, z 1 . 

For I « L, a r and A I * I - I „ we get for the principal axis torque 

y 7* X 

due to the gravity gradient in the inertially llxod O system (with 



cN 

0 

sn" 


0 


rg sN 

K K !• 

m gg: 1 

0 

1 

0 


~‘’3 

J5S K 1* 

gg f 

-*■3 


»sN 

0 

cN 



I'j ON 

M M 


Only the torque about the z a?ds is of interest 
^ggz “ K^^zsnsy(cycn-sysn) 

With n = Ks2y; sn ? Ks2yj on k I and K = -• we get 

" Kg^zKs2ysy (cy-syKs2y) 

a 0.25K^^zK [(1 - K) - (1 - K) c4y - 2yc4y) . 
EMmination of the cyclic terms loaves a constant bias torque of 


T_. a 0.25K ZK (1 ~ K) . 
HB 


Integrating over an orbit and substituting the value for 


z, K and K gives 

DO 


''HBB==Wns = - “o''' »yni "ym> 
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Inversion yields 


whore 

» 2H/[3fi^AI Hy„, (1 + Hy^)! , 

and e_„ is the normalized z momentum ramp» 
ssx 

Because of noise, a much smaller value had to be used in 
actual operation (see EOVV OPERATION and PERFORMANCE). 


APPENDIX E 


QUATERNIONS - A BRIEF EXPOSITION 


Complex numbers of the form z - a + ib have proven to be a valu- 
able concept in the study of many physical phenomena. A generalization 
of this concept which proves useful in the study of rotational motion is 
the quaternion. Recall that the imaginary unit i = Let us define 

additional units j and k together with the product operation o: 


ioisjoj=kok = “l 
ioi=-joi=k 


jok--koj=5i 

koi=“iok-j 


(El) 


We shall define a quaternion as any quantity of the form 


Q = Q4 4- i Q1 + j Q2 + k Q3 , (E2) 


(See Note 1 at end of Appendix E.) By analogy to the complex number 
terminology Q4 is referred to as the real or scalar part of Q. It will 
also be convenient to think of the remaining part of Q as the imaginary 
or vector part. The reason for this will become clear as we proceed. 
Let R = R4 + iRl + jR2 + kR3. The sum of quaternions Q and R is 
defined as 


S = Q + R = (Q4 + R4) + i (Q1 + Rl) 


+ i (Q2 + R2) + k (Q3 + R3) . (E3) 


From this definition we can see that the sum operation is commutative and 
associative. We can now give the complete definition of the product opera- 
tion o: 
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P ~ Q o R « (Q4 R4 » Q1 R1 - Q2 R2 » Q3 R3) 

+ 1 (Q4 R1 + Q1 R4 H- Q2 113 “ Q3 112) 

+ j (Q4 R2 Q2 R4 + Q3 Rl - Q1 R3) 

+ k (Q4 R3 + Q3 R4 + Q1 R2 ~ Q2 Rl) . (E4) 


With this dofinition we can show that o is associative and distributive but 
not commutative, i.o. , Q o R =it R o Q. Wo shall call any quaternion hav- 
ing zero imaginary part a scalar and, obviously, the algebra of scalars 
is just the algebra of real numbers. Thus, multipliealion of a quaternion 
by a scalar simply results in a quaternion whoso elements are multiplied 
by that scalar according to definition E4. We can now doJine the differ- 
ence operation as 


D ~ Q - R - Q -I- (“D 0 II . (E5) 


For convenience, we shall always omit the o when multiplying a quaternion 
by a scalar so that (-2) o Q - -2 Q. 

By analogy to complex algebra, lot us deiine the conjugate quater- 
nion to Q. The conjugation operation will be denoted by ( )*♦“. Thus 


Q* - Q4 - i Q1 “ j Q2 - k Q3 


(E6) 


So far all of our definitions have been extensions of those for complex 
numbers ns can be seen by assuming Q2 s Q3 s r 2 - E3 " 0. Thus the 
complex number system is a subset of the quaternions. It can be easily 
seen that 


Q* 0 Q ~ Q 0 Q* = Q1 Q1 + Q2 Q2 + Q3 Q3 + Q4 Q4 , (127) 


Note that Q* o Q is a pure real mimber or scalar. With tins observation 
we can define the inverse: 


^ (1/(Q* 0 Q)) Q*^/(Q* 0 Q) ; Q*o Q ^ 0 . (E8) 
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llnttlly then, wo can define a division operation: 


Q : R s Q o K 


(E9) 


We can sefe that Q : Q = Q o q"^ = o Q s 1 so thot Q**^ sotisfies the 
necessary properties of an inverse as long os Q 0. This will be useful 
lotcr. 

We need some additional results and definitions. First we can show 

that 


(Q o R)* = IX^ o q* 


(ElO) 


If the quaternion Q = Q*, then Q is necessarily a scalar. Also, if Q » 
“Q*, Q is purely imaginary or a vector quaternion. If V is a vector 
quaternion, wo shall designate this by an undoi’line as is also used to 
designate a 3-spaco vector, i.c., (i, j, and k will not be underlined) 


V = i VI + j V2 + k V3 


(Ell) 


For compactness of our notation we shall lot 


Q = Q4 t- Q 


(E12) 


where 


g = i Q1 + 3 Q2 + k Q3. 


Thus 


Q o R s (Q4 Q4 - g • R) + Q4 R I- R4 g + Q X R . (E13) 


The operations • and x ai’c defined as for 3-spaco vectors so that 


g • R = Q1 R1 + Q2 R2 + Q3 R3 


(E14) 


and 


g X R = i (Q2 R3 - Q3 R2) j (Q3 R1 - Q1 R3) 


+ k (Q1 R3 - Q3 .'ll) . (E15) 


For vectors A and B, 


AoB="A*B+AxB . (E16) 

In general the product of quaternions mixes scalar and vector parts 
together so that tliis product is not very interesting in the study of 
rotational motion in 3~apace. The triple product 


V = Q* o V 0 Q (E17) 


is more interesting since it does iJ\*eserve scalar and vector parts of V 
without mixing them. This property is trivial for the scalar part of V 
and follows for the vector part since 


V'* = (Q* o V o Q)* - -Q* o V 0 Q :s -V' 


(E18) 


Hence as noted the triple quateniion product (E17) takes a scalar into 
a scalar and a vector into a vector for any quaternion Q. Furthermore, 

the length of the vector |V | = /v • V and 


T • Y' = V' o y* = Q=»* o y o Q 0 (Q* o V o Q)* 


= (Q* o Q)^ y o y* 


(E19) 


Equation (E19) indicates the triple product (E17) multiplies vector length 
by the factor Q'^' o Q which is a real number. We note that if Q* o Q = 1, 
vector length is preserved and the vector mapping V ^ V’ looks like 
a rotation operator. It is a lineur operator in that aA +~bB *>■ aA' + bB'. 
Restricting ourselves to iiormalized quatei'nions wMcli' pi’esorve length, 
we sec that E17 is equivalent to a rotation of vector y into y. Since 
we are looldng only at normalized quaternions, we can without loss of 
generality, repi'esent Q as 
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(E20) 


Q = cos (|)/2 + sin ij)/2 u ; whore u • u ~ 1 


The triple product equation (E17) can be combined with equation (E20) 
to give 


V' = cos (}i V “ sin <l>u^V + (l - cos ‘f’) u u • V . (E21) 


Equation (E21) is the general form of the i-otation of a vector V about 
axis u througli the angle -0 . That this is true is soon by oxaniining the 
i-otation operation. Clearly any vector along the rotation axis u is not 
changed by the rotation so that if the vector V is broken into parts par- 
allel to and normal to u; i.e. “* 


V = V + V , whore V = V • u u and V • u ™ o . (E22) 

“ "’ll “ I "”11 “ “ “I * 

Wo must also have V' ~ + V'^ ; whore V' • u 0 and V « V' . 

Since V and V are normal to\i, wo can express V ns 
•X "X ‘‘ ~x 


V 

“i 


X u X V + y u X (vj X V ) 

" ““X ” " "X 


(E23) 


Now V 'V =1^1 eosa = ~y |V^ 1 y = -cos a and V x V = 

1 i^X. "X X “X “X ’"1 

|V^ I sin a u and thus x - sin a; where a is the rotation angle, Com- 
bining, wo obtain 


sill a u X V - cos a u x (u x V ) 

“1 - "X 


(E24) 


and 


“ H ■ Y E ^ Y “ * y. ^ ^ Y) • (lil25) 


The equivalence between equations (E21) and (E25) for a « - «j? is estab- 
lished. Thus, the mapping (E17) is equivalent to a rotation operator in 
a vector 3-space. Since the coordinate directions arc also vectors, we 
can ixitate the coordinate system instead of the vector. Rotating the 
vector through -ip yields tlie same components V' as rotating the coordi- 
nate axes through tj>. Thus we can look at V' as a new vector formed 
from V by rotation and expressed in the old lixis system or as the old 
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vector expressed in new axes rotated relative to the old. Wo have now 
demonstrated that any coordinate system rotation can be roprcsontod by 
a quaternion. Note that if Q satisfies equation (E17), then so does -Q. 
Looking back at equation (E20) tells us tliat -Q corresponds to (j> + 360 
degrees wJiich represents the same attitude ij» does. 

We now look at the time variations of Q. Since Q is constrdned to bo 
normalized, wo necessarily have 


d/dt(Q* o Q) - 0 s o Q + Q* o Q . (E26) 

We see from equation (E26) that Q* o “ (Q’*' o <^)* and hence must 
be a vector. Let us define tliis vector by 

Q* 0 Q s 1/2 u) (E27) 

-j- Q = 1/2 Q 0 w (since Q o Q* = Q* o Q ~ 1) . (E28) 

We shall see the reason for the 1/2 factor later. When we evaluate the 

rate of change of a vector with time in two reference frames, we find 

^ Q o V' o Q* + Q o V' o Q* + Q o V' o Q* . (E29) 

Using equations (E28) in (E29), 

V = Q 0 [V' + 1/2 (ui o V' - V' o u)3 o Q* (E30) 

y = Q 0 [V' + w X V'] o Q* . (E31) 


Now the I’eason for the factor 1/2 becomes clear. It is so that we can 
identify w. . Equation (E31) is exactly like the coiTesponding equation 
for 3-sp^e vectors if we identify « as the angular velocity of the primed 
reference frame expressed in prim^ coordinates. This identifcation fol- 
lows from the fact that eq\iation (E31) holds for an arbitrary vector V. 
Thus , w is identified as the relative angular velocity of the primed axes 
with respect to the unprimed. 


The above discussion completes the basic development of our qua- 
ternion tools. We now turn to the problem of developing a more conven- 
ient notation. The most logical choice which comes to mind is o matrix 
representation. The quaternion Q would logically become 


Q = 


Q3. 

Q2 

Q3 

Q4 


(E32) 


Looking back to the definition equation (E4) of the quaternion product 
o, we see that for P Q o R we have 



Q4 

-Q3 Q2 

Q1 


III 



Q3 

Q4 -Q1 

Q2 


R2 

"w 

p - 

-Q2 

QX Q4 

Q3 


113 

i\ Q R . (E33) 


-Q1 

-Q2 -Q3 

Q4 


R4 


(See Notes 

at end 

of Appendix E.) 

Similarly, i;. D = A o B o C, then 


D^AoBC=ABCsABC . (E34) 


This ri jult shows tlint the set of matrices of the form Q have the proper- 
ties. of the quaternions and in fact comprise a maxtrix repi’esentation of 
quaternion algebra with matrix multiplication corresponding to o. We can 
also express the quaternion product in the alternate form 


R4 R3 -R2 lU 


QX 

-R3 R4 RX R2 


Q2 

R2 -R1 R4 R3 


Q3 

-RX -R2 -R3 R4 


Q4 


(See Notes at end of Appendix E.) Here the mapping also yields a matrix 
representation except the order of the factors must be reversed. Thus, 
we have 


43 


P = Q R s 1 Q 

D«ABC-5 6a . (E36) 

wTW - c I 

According’ to these definitions and results, we have 

ABG-ACB“CBA-CAB . (1537) 

Thus, an intex’csting and sometimes useful result is that A U = C A, 

Let us now look at our previous work and make use of these new 
definitions : 


V' = Q* o V o Q ’>• y = Q* V Q ~ Q* Q V 


(1538) 


Note that we now liave a matrix formed form () which rotates coordinate 
axes and produces V' from V. Linear vector spaces are also represented 
by matrices. The 3»spaoo vectors, V and V arc related as 


V' “ M y 


(1539) 


Hero M is a 3 x 3 matrix which transforms components of V to primed 
coordinates. Referi’ing* back to equation (E25) and replacing n by 
we see that 


M 


u u 


sin i}> u cos (II - u u/') 


(E40) 


(See Notes at end of Appendix E.) The matrix u is the so~calIad cross 
prodiiet mati’ix and just liappons to be the upper 3x3 formed by dropping 

the final row and column of u. The matrix I is the identity matrix of the 
appi’opriate size to lit the current application. We have already shown 

the equivalence between M and Q* Q. Let us look more closely at the 
latter since it is 4 x 4. We can partition the double-tilde or double-bar 
mati’ices as 
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(E41) 


Q = Q4 a + 



isl 


r 1 **'1 

“Q 1 9 


; Q s Q4 11 + 

1 

( — 

l-9^' 

1 


T ' 

[-8 ! oj 


Since Q* o Q = 1, Q* = so that Q* Q = i. Also, Q’*' = Q (transpose) 
so that 


iw ^ 

Q* Q 


Q4^ 1 - 


2Q4Q'fQQ+QQ 

0 



(E42) 


Hence M = Q4^ 1 - 2 Q4 Q + 2 Q 1, In expanded foi’m 


9 9 

vjr * Q2“ 

i}3“ + 

1 

r j 2(01 (j2 + (j 3 on 

\ 2(Ql Q3 • Q2 04) 


2(Q2 Ql ^ 

i) iJ Q 4 ) 

j 01*'+ 02“ (}3“ + 04*^ 

1 

1 2(02 03 + Q1 Q4) 

(E43) 

2(Q3 i}J + 

Q2 

1 2(Q3 Q2 - Q1 Q4) 

i 

j Ql" • Q2^ + 03^ + Q4^ 



From vsquation (E28) 


(E44) 


Equations (E43) and (E44) summarize the useful results from our 
discussion. 

VVe are now ready to consider the question of successive I’otations 
applied to a cooi’dinate reference. A coordinate frame x’otation is a rigid 
displacement of all the poi* ts in the system with a fixed axis passing 
through the origin. Thus, it would seem that several successive rotations 
should displace every point except the oi’igin. Lot us now consider the 
coordinate frame as a rigid body and determine the most general displace- 
ment of it which keeps one point fixed. We must first explain what is 


Q 


Q o to — “g- (0 Q 


0 (03 -r»)2 wl 


Ql 

3 

o 

CO 


Q2 

(n2 -col 0 uC 


Q3 

”t«'l — ui2 — ui3 0 


Q4 
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meant by a rigid body displacement, A rigid body displacement is one 
which preserves distances between every possible poir of points in the 
body. The displacement is mathematically represented as a vector func- 
tion f. This function then has two basic properties: 


1) f(0) = 0 


2) |£(£a^ ’■ = IjLa " 


(E45) 


To study this in more detail, we define two additional points r^^ and rg 
together with their images under f, f(rp and fCrg). Let fj^ := f(rj^) and 
fg = define unit vectors 




j ~ 


£2 ~ £2 • II 
£2 - £2 • a 


- -2 ~ -2 * Hl-1 
l~2 “ -2 ' HiHl 


(E46) 


k = i X j ; Ug X Ug 


The vectors i, j, k and Uj^, Ug, Ug each form orthonormal bases for 3 
dimensional space. An arbitrary vector r can be expresses as 

r=xi + yj + zk . (E47) 

The corresponding f(r) = f u, + f„ Ua + f Uo. Condition 2 of equation 

X y z — u 

(E45) can only be satisfied if 

f(£) = X Uj^ + y Ug ± z Ug , (E48) 


Thus we can define two functions f+ and f- that both satisfy equation 
(E45) and .lap into f and rg into fg. The function f- can be viewed 
as the reflection (x,y,z) (x,y,-is) followed by f+. We are only 
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interested in continuous transitions fi*om an initial position to a final posi- 
tion and thus reflections must be eliminated since it is not possible to go 
fi*om (x,y,z) to (x,y,-z) continuously without violating condition 2 of 
equation (E45). Thus continuous rigid displacements can only occur in 
the form f+. This function can be written out ns 


f(r) = X u^ + y Ug + z: Ug = ct^i + ocyj + (E49) 

s X Uj^ * i y H2 ' ^ H3 * i 

Oy - X Uj^ • j + y Ug • j + 55 Ug • j (E50) 

-J ■ ^ y ^2 * ^ >*3 ' 


Equation (E50) can be rewritten in the matrix form 




(E51) 


The vectors a and r are of the same length and since this must hold for 
nil pairs a and r wq must have that 


- 1 


(E52) 


Equation (E52) also implies all eigenvalues of M are of unit magnitude. 
The eigenvalues and eigenvectors of M may be complex so that if M x = 
T* T T 

X X, then x M M x = 1 ~ X-i^Xx ’•‘x. Eor 3-dimensional space M must 
have at least one real eigenvalue. Since M is real, its eigenvalues must 
occur in complex pairs. Therefore at least one eigenvalue of M must be 
equal to 1. The value -1 could not bo acceptable since it would imply 
M u = “U whibh would be a reflection and already ruled out. Thus, we 
have that det M = 1. 

The matx'ix M is now looMng vei’y much like a rotation since the 
eigenvector u is an eigenaxis. All we must do now is to determine the 
angle of rotation. Along with the eigenvector u, let us deflne unit vec- 
toi’s y and w such that u, y, w is an orthonormal basis set. Also, we 
assume w ~ "u x y. Witir th’ese"" definitions we can express the matx'ix M 
as 




(E53) 
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From the fact that M u s u and that M = 1 equation (E53) reduces to 

T T T 9 9 

® H H. + P(l Z + iX \J!1 E ~ H X ) 5 p + q - 1 

(E54) 

We can now eliminate the vectors v and w from tliis equation by use of 
the proper function of u. Thus ”* ^ 


T T 

M = u u + p(l - u u ) - q u 


(E55) 


This completes the pi’oof that the matrix M is a rotation matrix. This is 
now obvious from inspection of equation (E55) by comparing it to equa- 
tion (E40) with p = cos Ip and q - sin if>. Thus the most general displace- 
ment of a rigid body (or transformation of a coordinate system) in wMch 
at least 1 point remains fixed Is n rotation about a fixed axis i,o., the 
final orientation can be obtained from the original by a single rotation 
about the axis u tlirough the angle p (u and p are determined from M) 
oven though the actual motion from initial to final may have boon more 
complex. 


What all the previous discussion boils down to is that the product 
of a pair of rotations is itself a rotation. Thus, if and Mg are rota- 


M, 


Mg is also a i*otation 


tions about u, and u, respectively, then M, 
through some angle p^ about some axis Ug, In fancier texmis the sot of 
rotations forms a group under matrix multiplication. 


The results of our previous discussions now suggest some new 
notation that may aid us in keeping up with the multiplicity of coordinate 
systems that must usually be doiUt with in analysis of spacecraft rotational 
dynamics. To remain completely general, let us consider tliree coordinate 
fx’amos A, B, C. We sh. ‘ let the symbol [BA] I’eprosent the rotation 
matrix which transfoxmis /ector expressed in tile A frame to a vector 
expressed in the B frame. 


= [BA] 


(E56) 


For convenience we use the notation superscript (A) or (B) ote. to indi- 
cate which coordinate frames the vectors ax’e being expressed in. If the 
superscripts are not specified, it means that the coox*dinate frame is 
impEcit in the definition of the symbol or that it doesn't matter as long 
as all vectors are in the same frame. We are here more interested in the 
rotations [BA] etc. Thex’e are three x’otations between pairs: AB, BC, 

CA. From our xxrovious work 
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[CAJ « tCB] [BA] 


(E§7) 


Corresponding to equation (E57) is a quaternion relation of similar form. 
First, since s; [CA] wo have 

~ ^CA* ® ^ ‘^CA ’ 

Hero, is the quaternion corresponding to [CA]. Thus analogous 
to equation (E57) wo have 

y«)=(Jj,,,*0Qy^*0V«)„Q^A0Qci, . (E50) 

Interestingly, wo see that o Q^jj so that the factors occur 

in reverse order from equation (E57). However, if wo use the double^ 
bar operator we can multiply in the same order, i.e. 


*^CA “ ^CB %A 


(EGO) 


This result is the one which we wish to use analogous to equation (E57). 
Equations (E57) and (EGO) have an easily remembered foi’ia and in fact 
behave as if multiplication cancelled the terms appearing on the inside. 
This makes it quite easy to construct chains of transformations to any 
desired system. lii this notation we see that 


[BA] = [AB]'**^ s; [AB]“^ ; also 


*^BA " ^AB* " ^AB 


(EGl) 


Finally, thex’C are some useful tidcks with the new notation we have 
defined, Referidng to equation (E28) and adding the subscripts we have 

defined, we have s 0.5 o The vector is the 

angular velocity of B relative to A witli components in B. Consider the 
quatei’nion Qgg. 
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^OB * I ^CB ° ® I ^CB j 

“ I ^CB ° ” I ^ ^‘^CB* ° i^BA^^^ ® ^CB^ 

-In ^ (C) 1 (B) ^ 

“ f ^CB ® i^’CA “ 2 “BA ° ‘^CB 

or in matrix form 

^CB " I ^CB ’ 

The utility of oqiiation (E62) is most apparent when wo use it to compute 
the attitude error of a spucocraft relative to a moving or moveable refer- 

once. Note that the components of anB Qi’o expressed 

in different frames. Normally, would come from rate sensors 

which are body fixed wliile is a commanded maneuver rate which 

is naturally defined in the moveable reference. Tins equation then allows 
us to use both quantities directly without either being transformed. 

It often becomes necessary to compute the quaternion corresponding 
to a given rotation matrix, i.e., find Q given [BA]. We have developed 
a computer algoi’ithm to do this. 


1) Define matrix 



A is the given rotatic i. 


’ 2) S' - S + + (1 - tr A) 1. 

I 
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3) I " max S'y (index of largest element along diagonal of S’). 

4) Qj = S'j^/2 /s^ . 

1 ~1 for Q4 < 0 

. 

+1 for Q4 > 0 


Another useful and perhaps obvious technique is the expression of the 
quaternion resulting from a sequence of Euler rotations (rotations about 
coordinate axes) : 



where s A sin and c A cos. The corresponding rotation is tBA] and is 
given by 


[BA] 


[tl* ]• 


n 




[<{> 


I'i. 


The vectors u can bo any of the three coordinate axes [1,0,0] , [0,1,0] 
or [0,0,1]^. If u= [1,0,0] then i^^ = 1, etc. We have added the con- 
vention that a x’otation braclcet with a subscript is an Euler rotation about 
the indicated axis. As an example consider the quaternion formed when 
ij^ — 1, i2 ~ 2, ig ^ 3; 
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■0 


‘<’1 
s 


c —g* 




+ G “2* 


<}»2 

s 


^^3 

T 


) 


+ j 


(• 


‘J’l 

T 


cl>2 

T 


t|.3 

T 


- s 


h 

T 


•^2 


) 


+ k 


( 


. *1 
c 


« ‘^’2 
c " 2 " 


2 


+ s 


^2 


T^) 


the cowesponding [BA] is 


1 

C().2Ci|)3 { 

stj) + ci|)j^S(j)3 1 

- Cc|l j^SfJlgCiltg + 


1 

- C+jStj 1 

- sij)j^s4>2«‘l’3 + ct|ij^C(j)3 j 

c«J>^st|i2Si|>3 + 

Si|»j^C(t>3 

^■l>2 i 
1 

- s4,j04,2 1 

cit!j^ci|>2 



In this brief exposition, we have developed a number of useful quaternion 
results and notations. This by no means exliausts the possibilities. The 
available quaternion litcratui’c does not present the material in an easily 
applicable form and thus this short development is presented to fill that 
g«P- 


NOTES: The following notes apply to the previous discussion: 

1) We use Q4 rather than QO for convonience. Since these quater- 
nion equations will be adapted foi' the computer and since 0 is not usually 
allowed as a subscript it becomes necessary to use something else. We 
desire to use 1, 2, 3 for the vector components, hence Q4 is the real 
part; 

2) The symbol ~ is called double tilde and the symbol =^^3 is called 
double bar; 

3) We shall define 



-Q3 

0 

Q1 


Q2 

-Q1 

0 


wliich is the tilde or cross product matrix for the 3-vector Q. 
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APPENDIX F. ROLL COMMAND 

X>I I 


To uso tho acquisition sun sonsoi* (ACQ SS) for a two axes strapdown 
update, tho sun Uno has to nominally pass through tho contor ol' tho 
ACQ SS onco per orbit, no matter what tho tilting angle is. This can bo 
done by rolling tho vehicle about tho pi^ncipal x axis by tho angle n 

AkI I 


such that the veMcle z axis always has an elevation angle ii above the 

iV 

orbital plane, equal to the elevation angle ij of tho sun. Tho transfor- 

A 


mation 

from 

tho 

L system 

to. the V system is 

(s " sin, 

e « cos) 



'Ml 

'''ii 

'Mu’ 


1 

0 0 


.11 



VII 

0 

• 

” yii 

[VI. 1 

'Ml' 

'M';i 



I) 

' vi» '* \n 


.11 

‘•..11 " 


0 

1 

0 






II 

Ml Ml 


ll 

ll ! 


Mvii " 

' VII. 


from which wo derive 


VL 


32 


sn,,„ + 1 C 23 


K33 sn_^„ 0,1. 


‘zn 


or 




I 


with 
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This yields 


T] 


XU 


arc tan 


jKog/K jjjj I 


•r aro! tan 


n //I 
'xp' 


■xp 
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Inverse tangents ore used since no other ARC functions are available on 
bo»u*d. The values of the K's ore given in APPENDIX G. They apply 
for the EOW A orientation. For the EOVV B orientation wo have 


[K] 


i.o.. Kj3 


B 


•10 0 
0-10 
0 0 1 


tK] 


and Kgg change sign. 


APPENDIX G. DATA 


Vehicle nwments of inertia matrix 

(kgm^) 


894828 

-63414 

-529360 


-63414 

3763111 

-27295 


“529360 

“27295 

3598005 


Principal moments-of-inertia matrix (kgm^) 


793332 

0 

0 

0 

3767879 

0 

0 

0 

3694732 


Transformation from veliicle to principal coordinate system 


[ig = [PV] == 


0.982357 

0.017288 

-0.186213 


0.022682 0.185632 
0.977351 -0,210913 
0.210401 0.959716 


















LO: 


»Jin “ - * db) 

^K13* ( IF EOVVBFLAC3-0 

I -0,1tS633 IF EOVV B FLAG ¥= 0 

tj*p - (tin Tj^ + Eki 3 iJn / (0.BB2619 cm »?,„) 

riun " ARCTAN <r?^p / (1 - 

t?x„ + 12.396® IF EOW B FLAG =A 0 


(in 


iln’'*"/2 

,co(’^*"/2 


O 

(in ’'V"/2 


, ’?yn/2 


°BL * Qm Qn 

M “ 2 IqBn qgL2 + 9BL3 RBLd^ 

'2 " <‘>iL2 «}iL1 ~ qiL3 + ‘»iL4> 

f3 * 2(qBL3 'IBL2 ~ tlBLI ‘»BL4> 
nb*0~2)?y^ co( izn^ t^^)) 

®x - rifio^b 

© y = f2 ^^b 

• 

® z * ''3 ^b 


°BL QlI 



STRAPDOWN UPDATE REQUIREMENTS 




RSSP •» RSS 


RSS < B.0“ 
\ 7 ✓ 


RSS RSSP 


« 

“COI Jm ygy 

$2 ® 

•in 

$3 o 

COI coi-j'gy 

V, . 

2 (qvii tlvi3 ”RV(2‘«VI4l 

V 2 - 

2 RVl-J +‘»V(2‘>VI3> 

V 3 = 

<9VI^3 -‘Iv |^2 

AQl = 

0.5 (-S 3 V 2 t SjVa + ^^s,) 

AQj = 

0.5 (S 3 V, "SiV 3 + Ai^S 2 ) 

AQ 3 » 

0.5 (-S 2 V, +S,V 2 + M,S 3 ) 

tl 

0 

< 

(1 -AQ,2 -AQ22-AQ32)!^ 

(ovi) 

•> ( AQ ) { Qv, ) 

* 0 ; 

^yc “ ’^yc + Trij,* ^ ^2 ^°2 *' ^3 AQ 3 ) 
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